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Abstract 

We study prompt J ftp— meson hadroproduction invoking the hypothesis of gluon Reggeization 
D . in t— channel exchanges at high energy and the factorization formalism of nonrelativistic quantum 



chromodynamics at the leading order in the strong-coupling constant a s and the relative velocity 
of quarks v. The transverse-momentum distribution of direct and prompt J/ip— meson produc- 
tion measured at the Fermilab Tevatron fitted to obtain the nonperturbative long-distance matrix 
elements, which are used to predict prompt J/ip production spectra at the CERN LHC. At the 
numerical calculation, we adopt the Kimber-Martin-Ryskin and Blumlein prescriptions to derive 
unintegrated gluon distribution function of the proton from their collinear counterpart, for which 
we use the Martin-Roberts-Stirling-Thorne set. Without adjusting any free parameters, we find 
good agreement with measurements by the ATLAS, CMS and LHCb Collaborations at the LHC 
at the hadronic cm. energy y/S = 7 TeV. 
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I. INTRODUCTION 



The production of heavy quarkonium at hadron colliders provides useful laboratory for 
testing the high-energy limit of quantum chromo dynamics (QCD) as well as the interplay 
of perturbative and nonperturbative phenomena in QCD. 

The total collision energies, \/S = 1.8 TeV and 1.96 TeV in Tevatron Runs I and II, 
respectively, and v^S* = 7 TeV or 14 TeV at the LHC, sufficiently exceed the characteristic 
scale [i of the relevant hard processes, which is of order of quarkonium transverse mass 



M T = y M 2 + p?, i.e. we have Aqcd <C /i <C y/S. In this high-energy regime, so called 
" Regge limit" , the contribution of partonic subprocesses involving t-channel parton (gluon or 
quark) exchanges to the production cross section can become dominant. Then the transverse 
momenta of the incoming partons and their off-shell properties can no longer be neglected, 
and we deal with "Reggeized" t-channel partons. These t— channel exchanges obey multi- 
Regge kinematics (MRK), when the particles produced in the collision are strongly separated 
in rapidity. If the same situation is realized with groups of particles, then quasi-multi- Regge 
kinematics (QMRK) is at work. In the case of J/ip— meson inclusive production, this means 
the following: J/ip— meson (MRK) or J/^-meson plus gluon jet (QMRK) is produced in the 
central region of rapidity, while other particles are produced with large modula of rapidities. 

The parton Reggeization approach (PRA) is particularly appropriate for high-energy 
phenomenology. We see, the assumption of a dominant role of MRK or QMRK production 
mechanisms at high energy works well. PRA is based on an effective quantum field theory 
implemented with the non-Abelian gauge-invariant action including fields of Reggeized glu- 
ons 3] and Reggeized quarks jjj. Reggeized partons interact with quarks and Yang-Mills 
gluons in a specific way. Recently, in Ref. j^J, the Feynman rules for the effective theory of 
Reggeized gluons were derived for the induced and some important effective vertices. This 
approach was successfully applied to interpret the production of isolated jets [6j, prompt 
photons 7], diphotons 8], charmed mesons js], bottom-flavored jets 10] measured at the 
Fermilab Tevatron, at the DESY HERA and at the CERN LHC, in the small-pr regime, 
where pr « \/S. We suggest the MRK (QMRK) production mechanism to be the dom- 
inant one at small pt values. Using the Feynman rules for the effective theory, we can 
construct heavy quarkonium production amplitudes in framework of non-relativistic QCD 
(NRQCD)[u, 12]. 
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The factorization formalism of the NRQCD is a rigorous theoretical framework for 
the description of heavy-quarkonium production and decay. The factorization hypothe- 
sis of NRQCD assumes the separation of the effects of long and short distances in heavy- 
quarkonium production. NRQCD is organized as a perturbative expansion in two small 
parameters, the strong-coupling constant a s and the relative velocity v of the heavy quarks 
inside a heavy quarkonium. 

Our previous analysis of charmonium jl^, and bottomonium fl^ . \i\ production at 



the Fermilab Tevatron using the high-energy factorization scheme and NRQCD approach 
has shown the efficiency of such type of high-energy phenomenology. In this paper we repeat 
our calculations for the direct and prompt J/ip— meson transverse momentum spectra at the 



Fermilab Tevatron 



13[ to obtain by fitting procedure octet nonperturbative matrix elements 



(NMEs), than we calculate prompt J/ip— meson spectra, which were measured recently 
at the CERN LHC Collider at the energy of y/S = 7 TeV. We find a good agreement 
of our calculations and experimental data from ATLAS 16], CMS 17] and LHCb [lsl ] 
Collaborations. 



II. MODEL 

Working at the leading order (LO) in a s and v we consider the following partonic sub- 
processes, which describe charmonium production at high energy: 

R( qi ) + R(q 2 ) -+ hM\*S?\ 1 S®,>P?]( P ), (1) 
R( qi )+R(q 2 ) -+ HFS?](p) + gtf), (2) 

where R is the Reggeeized gluon and g is the Yang-Mills gluon, respectively, with four- 
momenta indicated in parentheses, %\n\ is the physical charmonium state, n = 2S + 1 £ J y 8 ' 
is the included cc Fock state with the spin S, total angular momentum J, orbital angular 
momentum L and in the singlet ^ or in the octet ® color state. 

In the general case, the partonic cross section of charmonium production receives from 
the cc Fock state [n] = [ 2S+1 L ( j' 8 ' ) ] the contribution 
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d Z(R + C c[^L { j l > 8) ] ^U) = da(R + R-> cc[ 2S+1 L^}) i ^ 1 ^]) ; (3) 

A C olA pol 
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where iV co i = 2N C for the color-singlet state, iV co i = N 2 — 1 for the color-octet state, and 
A^ po i = 2J + 1, (O u [ 2S+1 L { j' 8) ]) are the NMEs. They satisfy the multiplicity relations 

( lKnS)[3p(8)]) = (2J+ l)(O^" 5 )[ 3 P (8) ]), 

(O x ^Pf ] \) = (2J+1)(0**[ 3 P O (1) ]), 

(O x ^[ 3 S{ 8) )) = (2J+ 1){0™[ 3 S[ 8) )), (4) 

which follow from heavy-quark spin symmetry to LO in v. 
The partonic cross section of cc production is defined as 

da(R + R -> cc[ 2S+1 4 1,8) D = + ->• ccp^L^DI 2 ^, (5) 

where / = 2x\x 2 S is the flux factor of the incoming particles, which is taken as in the 
collinear parton model 19], A(R + R — > cc[ 2S+1 l[j ,s ^]) is the production amplitude, the 
bar indicates average (summation) over initial-state (final-state) spins and colors, and rf$ 
is the invariant phase space volume of the outgoing particles. This convention implies that 
the cross section in the high-energy factorization scheme is normalized approximately to the 
cross section for on-shell gluons in the collinear parton model when q 1T = q$T = 0. 

The LO results for the squared amplitudes of subprocesses ([TJ and (J2J) that we found 
by using the Feynman rules of Ref. ||j] coincide with those we obtained in Ref. The 
formulas for the squared amplitudes \ A(R + R — > cc[ 2S+1 Lj 1 ' 8 '*])] 2 for the 2 — > 1 subprocesses 
([I]) are listed in Eq. (27) of Ref. Jjij]. The analytical result in case of the 2—^2 subprocess 



J2]) is presented in Ref. 



14j, where the results for the 2 — > 1 subprocesses are also listed, but 



in another equivalent form. The relation between these forms is discussed in Ref. 151 ] . 

Exploiting the hypothesis of high-energy factorization, we may write the hadronic cross 
section da as convolution of partonic cross section da with unintegrated parton distribution 
functions (PDFs) <3>^(ic,t, /i 2 ) of Reggeized gluon in the proton, as 

da( P + P ^n + X) = J d f f*^ Vg ( Sl ,t ul t) 

^ g {x 2 ,t 2 ,fj, 2 )da(R + R^n + X). (6) 



x 



^i = |qir| 2 ; t 2 = | q2T 1 2 ? %i an d x 2 are the fractions of the proton momenta passed on to the 
Reggeized gluons, and the factorization scale fi is chosen to be of order Mt- The collinear 
and unintegrated gluon distribution functions are formally related as 

xG p (x,fx 2 )= $ p Jx,t,/j?)dt, (7) 



so that, for qi^ = q 2 T = 0, we recover the conventional factorization formula of the collinear 
parton model, 

da(p + p^H +X) = J dx^ix^ /i 2 ) J dx 2 G p (x 2 , fi 2 )da(g + g^U + X). (8) 

We now describe how to evaluate the differential hadronic cross section from Eq. §5§ 
combined with the squared amplitudes of the 2 — > 1 and 2 — )• 2 subprocesses ([1]) and (T5]), 
respectively. The rapidity and pseudorapidity of a charmonium state with four-momentum 

_ p T; p 3 ) are given by 

1 p°+p 3 1 |p| + p 3 

y = n ln ~o 3' ^ = o ln T~\ 3' 9 

2 p u — p 6 2 |p| — p 6 

respectively, and dy = ^-drj. 

The invariant phase volume c?$ in the Eq. (E]) for 2 — > 1 subprocess ([1]) can be presented 
as follows: 

d 3 p 



d$(p) = (27r)V 4 )(g 1 + g 2 -p) 



(27r) 3 2p ( 



o 



^PTxfC Z^+^Wt P°~P\x2f , \j J /mt 

= — — 5 te 5 fe j=-)6 (qir + Q2T - p T )dp T dy. (10) 

From the Eqs. (jSJ), © and ffTUj) we obtain the master formula for the 2 — > 1 subprocess ([T]): 



dt\ J dip i 



dp T dy {p\ + M 2 ) 2 

x $^(6, t ls //)<^(6, t 2 , /i 2 ) ^(i? + R ->■ | 2 , (11) 

where t 2 = £i + p 2 - — 2pr-y/ii cos(0i) and the relation ^ 2 S = Pt + ^ 2 nas been taken into 
account. 

The invariant phase volume <i$ in the Eq. (JSJ) for 2 — >• 2 subprocess © can be presented 
as follows: 

d 3 p d 3 p' 



d$(p,p') = (27r)V 4 )(g 1 + g 2 -j9-p / ) 



(27r) 3 2p° (27r) 3 2p'° 



^((g 1 + g 2 -p) 2 )^ r ^. (12) 



Such a way, accordingly the Eqs. (JSJ), (JBJ) and f|T2|) . we have the master formula for the 
2 — ^ 2 subprocess (j2J): 

da(p + p -»■ H + X) p T 



y cz^i y c?x 2 y dt 2 y g^ 2 



^^'''^^''iS^- « i3) 



where 01,2 are the angles enclosed between qi^T and the transverse momentum pV of H 

1 



Xl 



,qiT + Q2T — Pt) 



M 



(14) 



In our numerical analysis, we adopt as our default the prescription proposed by Kimber, 
Martin, and Ryskin (KMR) 20[ to obtain unintegrated gluon PDF of the proton from the 



conventional integrated one, as implemented in Watt's code 2l|. As is well known 22], 
other popular prescriptions, such as those by Blumlein 23| or by Jung and Salam 24j . 
produce unintegrated PDFs with distinctly different t dependences. In order assess the 
resulting theoretical uncertainty, we also evaluate the unintegrated gluon PDF using the 
Blumlein approach, which resums small- a; effects according to the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) equation [sj]. As input for these procedures, we use the LO set of the 
Martin- Roberts-Stirling-Thorne (MRST) 25|, |26| proton PDF as our default. 

Throughout our analysis the renormalization and factorization scales are identified and 
chosen to be \x = £Mt, where £ is varied between 1/2 and 2 about its default value 1 to 
estimate the theoretical uncertainty due to the freedom in the choice of scales. The resulting 
errors are indicated as shaded bands in the figures. 



III. RESULTS 



At first, we examine direct and prompt (sum of direct production, production via radiative 
decays of Xcj mesons and production via decays of ip' mesons contributions) J/ip— meson 



production in proton-antiproton collisions at the Tevatron at the energies yS = 1.8 TeV 



(Run 1)0 and y/S = 1.96 TeV (Run II)(28j in the central region of pseudorapidity \r]\ < 0.6. 
The data of Run II includes pt distributions of prompt J/ip mesons, so far without separation 
into direct, Xcj— decay, and if)'— decay contributions. In the Ref. 13|, we have performed a 



joint fit to the Run I and Run II CDF data 



27 



28 



to obtain the color-octet NME's for 



J/ip, ip', and Xcj mesons. Our fits included five experimental data sets, which come as pt 
distributions of J / ip mesons from direct production, prompt production, \ r j d ecays, and tp' 
decays. Since in the previous calculations we have used the old set MRST98 25[ as a collinear 
input for unintegrated PDF, in the present study we repeat our fit with the next-generation 
MRST set 26[. We find small differences between the old and new fit parameters, however, 



it is important for precise description of the data. The results of our fit are presented in 
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the Table I along with results of the fit in the next to leading order (NLO) of collinear 
parton model (PM) and NRQCD approach 31] . Oppositely the Ref. 3l|, we perform a fit 
procedure by assumption for NMEs to be only positive. Than, using the CDF data for 
a prompt J/if> production 27], presented separately for direct J/ip mesons, J/if> from if)' 
decays, and J/ijj from Xcj decays, we obtain color-octet NMEs (O j ^^S[ 8 \ 1 Sq , 3 Pq 8 ^]), 
(0^'[ 3 s[ 8 \ l Sf 3 P^ 8) }), and (0^[ 3 S[ 8) }) independently from each other. 

Looking at the Table I, we find a good agreement with the NLO fit in collinear parton 
model performed in the Ref. 3l|, which strongly improves if we take into account that a 
sum of contributions of NMEs (O^^S^]) and (0 J ^[ 3 P O (8) ]) from the Ref. |3l|, leading to 
almost parallel J/ip transverse momenta spectra, corresponds to our contribution of the NME 
(O^pS?]). Such an agreement demonstrates a validity of factorization in the charmonium 
production in hadronic collisions, i.e. an independence of the cc production mechanism from 
the nonperturbative charmonium formation at the last step. It is necessary to note that a 
same consent between LO results obtained in the uncollinear factorization scheme and NLO 
results obtained in the collinear parton model is also observed when describing other relevant 

nn 

processes, see Refs. |6H10]. 

In Figs. 1-4, we compare the CDF data on J j ip mesons from direct production, ip' decays 



n 



28 



with the respective 



and Xcj decays in Run I [27] and from prompt production in Run II 
theoretical results evaluated with the NMEs listed in Table I. As default, we present in all 
figures the theoretical results which are obtained using KMR unintegrated gluon density 



20] . The comparison between KMR [20] and Blumlein 23] PDFs is made in Fig. 4 for a 



prompt J /if) production only. 

In Fig. 1 one can find a dominance of color-octet contributions at all values of direct J / if) 
meson transverse momentum: [ 3 5i ] contribution dominates at the large values px > 10 
GeV, and [ 1 5'q 8 ' ) ] — at the small px < 10 GeV. The situation is very similar for J/if) produc- 
tion from if>' decay, considered in Fig. 2. It is also important, that the [ 3 -Pj ] contribution 
vanishes in direct J / if) and if)' production. The obtained results are in agreement with previ- 



ous calculations of Ref. 



13] with a slight difference. In Ref. 29] it was also shown that in the 



direct J/if> production at the Tevatron color-octet contribution dominates. Oppositely our 



conclusions, in Ref. 



the main contribution comes from the [ 1 5'q 8 ' ) ] NME at all transverse 



momenta. The reason of this discrepancy arises from the fact that in Ref. 



29] the color-octet 



NMEs for J j if) meson have been obtained by a fit of direct J / if) production data for px > 5 
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GeV only. In our fit we take into account both direct [27| and prompt [28| J /ip production 
data, the last ones contain points in a small pr region. We observe, the inclusion of prompt 
J/ip production data in the fit to change the relative weight of different color-octet NMEs. 
This fact should be important in study of polarized J/ip mesons. 

In case of Jjip production via decays of Xcj mesons, considered in Fig. 3, we confirm 



the conclusion of Refs. 



13 



30j , which reads, that in the high-energy factorization scheme, 



the color-singlet contribution is sufficient to describe the data for production of P— wave 
charmonium states. 

In Fig. 4, the pt distribution of prompt J / ip production in Tevatron Run II is presented 
as a sum of the contributions from direct production, ip' decays, and Xcj decays. We observe 
at the pt > 2 GeV the J/ip mesons to be produced preferably directly. The contribution 
from Xcj decays dominates at small px < 2 GeV. The contribution from ip' is smaller than 
other ones and it exceeds the contribution from Xcj decays only at the pt > 16 GeV. 



The curve number (5) in Fig. 4 is obtained using Blumlein unPDF i 



The visible 



difference between this curve and the curve (4), which is obtained using KMR unPDF [20!], 
takes place only in the region of small J/ip transverse momentum. In the range of Pt > 5 
GeV, there is no difference between them as in the prompt production as in the direct 
production or in the production via decays of high charmonium states, ip' and Xcj- 

As it is obvious from Figs. 1-4, the theoretical uncertainties associated with the variation 
of the factorization scale \i are large at the small px region, taking a value of about factor 
5 between upper and lower boundaries, and they sufficiently decrease down to a factor 2 at 
the px > 6 GeV. The uncertainties from errors in the color-octet NMEs are small, they are 
about 7-10%. 

Moving on from Tevatron to the LHC, which is currently running at the total energy being 
about 3.5 times larger than at the Tevatron, we expect the range of validity of our approach 
to be extended by the same factor, to pt < 70 GeV, as we describe well the Tevatron data 
at the range of < px < 20 GeV. This expectation is nicely confirmed in Figs. 5-6, where 
the recent measurements of the prompt J jip production by the ATLAS Collaboration at the 



CERN LHC 



16], which cover the kinematic region 1 GeV < pt < 70 GeV and \y\ < 2.4, are 



compared with our predictions based on the particle Reggeization approach and NRQCD 



formalism. The measurements of the CMS Collaboration 



171 ] were performed in the similar 



kinematic range 6.5 GeV < pt < 30 GeV and \y\ < 2.4, see Fig. 7. We observe a dominant 
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role of direct production mechanism in the prompt J/ip hadroproduction at the all values 
of J/ip meson transverse momentum. Concerning the relative contributions of ip' decays 
and Xcj decays into a prompt J/ip production, we found the contribution from ip' decays to 
dominate at the large px > 20 GeV, and the contribution from Xcj decays to dominate at the 
small pt, respectively. Additionally, we compare our predictions with the data from LHCb 
Collaboration 18|, which were extracted in the range < pr < 14 GeV and 2 < \y\ < 4.5. 
We find a good agreement between our predictions and prompt J / ip production data at the 
moderate rapidity interval 2.0 < \y\ < 3.5, see Figs. 8-9. At the same time our theoretical 
result overestimates the data of at most factor 2 in the range of large rapidity 3.5 < \y\ < 4.5. 
This distinction is expected in the parton Reggeization approach, because the multi-Regge 
kinematics conditions to be broken if J/i[) mesons are produced with large rapidity. 

We observe, in Fig. 10, that relative contributions of the color-singlet (curve 1) and color- 
octet (curve 2) production mechanisms to the prompt J / ip spectrum strongly depend on the 
J ftp transverse momentum. Similarly to the NLO calculations in the collinear parton model, 
the color-octet contribution dominates at the large pt region, basically via the contributions 



of the co 
data 



16 



ior-octet NMEs (O j ^[ 3 S[ 8) \) and (0^'[ 3 S[ S} }}. It is significant, the experimental 



3 C (8)l 



- |l8| depend on the assumption of polarization of produced J/ip mesons slightly. 



We perform calculations and make a comparison to the data in a case of non-polarized J / ip 
meson production. 

Comparing our results with the recent studies of J/ip meson hadroproduction in the 
conventional collinear PM, which were performed in full NLO approximation of NRQCD 
formalism 3l| or in the non-complete NNLO* approximation of color-singlet model 32j, we 
should emphasize the following. At first, oppositely to NLO and NNLO* calculations, which 
provide a good description of data only at non-small pt > 5 GeV, we can reproduce data 



well at all transverse momenta p?- At second, the present study 
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13 



15 



along with the previous 



29 



30] demonstrate the 



investigations in the parton Reggeization approach [6- 
important role of (quasi)multi-Regge kinematics in particle production at high energies, 
this feature is out of account in the collinear PM. Such a way, we find the approach based 



and high-energy factorization scheme 



on the effective theory of Reggeized partons 0, Q 
with unintegrated PDFs, which in the large logarithmic terms (ln(/i 2 /Aq CD ), ln(Sy/^ 2 )) 
are resummed in all orders of strong coupling constant a s , to be more adequate for the 
description of experimental data than fixed order in a s calculations in the frameworks of 



9 



collinear PM. 



IV. CONCLUSIONS 



The Fermilab Tevatron and, even more so, the CERN LHC are currently probing particle 
physics at terascale cm. energies \/S, so that the hierarchy Aqcd <C /•* ^ which 
defines the MRK and QMRK regimes, is satisfied for processes of heavy quark and heavy 
quarkonium production in the central region of rapidity, where /i is of order of their transverse 
mass. In this paper, we studied QCD processes of particular interest, namely prompt J/ip 
hadroproduction, at LOs in the parton Reggeization approach and NRQCD approach, in 
which they are mediated by 2 — >■ 1 and 2 — > 2 partonic subprocesses initiated by Reggeized 
gluon collisions. 

We found by the fit of Tevatron data that numerical values of the color-octet NMEs are 
very similar to ones obtained in the full NLO calculations based on NRQCD approach 31 |. 
Using these NMEs, we nicely described recent LHC data for prompt J/ip meson production 
measured by ATLAS 16], CMS 17] and LHCb [l8] Collaborations at the whole presented 
range of J j ijj transverse momenta. We found only one exclusion, the region of large modulo of 
rapidity \y\ > 3.5, where LHCb data are by a factor 2 smaller than our predictions. However, 



this kinematica! 



and in Refs 



.0 



re gion is out of the applicability limits of the MRK or QMRK pictures. Here 



10 



13 



15 



29 



30J, parton Reggeization approach was demonstrated to be a 



powerful tool for the theoretical description of QCD processes in the high-energy limit. 
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TABLE I: NMEs for Jhb 



colline ar p artem model 
KMR 



tp', and Xcj mesons from fits of the CDF data 27J, [28J] in the NLO 



23j, and 



3l[ and in the parton Reggeization approach using the Bliimlein 
20l | unintegrated gluon distribution functions. The errors on our fit results are determined 
by the varying in turn each NME up and down about its central value until the value of x 2 is 
increased by unity keeping all other NMEs fixed at their central values^ When we obtained the 
value of x 2 /d.o.f. > 1, we have used normalizing multiplier approach 
unity. 



331 ] to reduce this value to 



NME 



PM NLO [21] 



Fit B 



Fit KMR 



(O^fsfty/GeV 3 
(0 J ^[ 3 s[ 8) ])/GeV 3 
(O^S^/GeV 3 
(0 J A/fP o (8) ])/GeV 5 



1.3 

(1.68 ±0.46) x 10~ 3 
(3.04 ±0.35) x 10~ 2 
(-9.08 ± 1.61) x 10" 3 



1.3 

(1.89 ±0.27) x 10~ 3 
(1.80 ±0.25) x 10~ 2 




1.3 

(2.23 ±0.27) x 10~ 3 
(1.84 ±0.19) x 10~ 2 




X 2 /d.o.f 



1.0 



1.0 

6.5 x 10" 1 
(9.33 ± 1.62) x 10" 4 
(3.27 ± 1.44) x 10~ 3 




(0^[ 3 s{ 1] ])/GeV 3 
(0^'[ 3 5j 8) ])/GeV 3 
(C^'[ 1 5S 8) ])/GeV 3 

(O^'[3p ( 8 )])/ GeV 5 



6.5 x 10" 1 
(1.88 ±0.62) x 10~ 3 
(7.01 ± 4.75) x 10~ 3 



6.5 x 10" 1 
(6.72 ± 1.15) x 10~ 4 
(3.63 ± 1.40) x 10~ 3 




X 2 /d.o.f 




0.033 


0.051 


[ 3 P (1) ])/GeV 5 
(O^ [ 3 5j 8) ])/GeV 3 


8.9 x 10~ 2 


8.9 x 10~ 2 
(2.14 ±0.67) x 10~ 4 


8.9 x 10~ 2 
(1.69 ±0.9) x 10~ 4 


X 2 /d.o.f 




0.89 


0.41 
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I — I 

FIG. 1: Direct J/ip transverse momentum spectrum from CDF Collaboration [27(, \/S =1.8 TeV, 
\q\ < 0.6, (1) is the contribution of [ 3 s[ 8) ] state, (2) - [ l S^ 8) ], (3) - [ 3 s[ 1] ], (4) - sum of their all. 
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FIG. 2: J/tjj transverse momentum spectrum from tp' decays from CDF Collaboration [27|, y/S = 
1.8 TeV, I17I < 0.6, (1) is the contribution of [ 3 s[ 8) ] state, (2) - [ 3 s( 1} ], (3) - [ l S^\ (4) - sum of 
their all. 



1(3 




FIG. 3: J/ip transverse momentum spectrum from Xcj decays from CDF Collaboration 



1.8 TeV, \r]\ < 0.6, (1) is the contribution of [ 3 P (1) ] state, (2) - [ 3 PH, (3) - [ 3 P 2 W ], (4) 



r3p(i)i 



N 8) ] 



(5) - sum of their all. 
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FIG. 4: Prompt J/tp transverse momentum spectrum from CDF Collaboration [28], y/S = 1.96 
TeV, \y\ < 0.6, (1) is the direct production, (2) - from Xcj decays, (3) - from tp' decays, (4) - sum 
of their all (KMR unPDF), (5) - sum of all contributions (Bltimlein unPDF). 
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FIG. 5: Prompt J/ip transverse momentum spectrum from ATLAS Collaboration [16], yS = 7 
TeV, (1) is the direct production, (2) - from Xcj decays, (3) - from ip' decays, (4) - sum of their 
all. For the different range in the rapidity: (a)- \y\ < 0.75, (b) - 0.75 < \y\ < 1.5. 
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FIG. 6: Prompt J/ip transverse momentum spectrum from ATLAS Collaboration [16], y/S = 7 
TeV, (1) is the direct production, (2) - from Xcj decays, (3) - from ip' decays, (4) - sum of their 
all. For the different range in the rapidity: a) - 1.5 < |y| < 2.0, (b) - 2.0 < |y| < 2.4 . 
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FIG. 7: Prompt J/ip transverse momentum spectrum from CMS Collaboration , VS = 7 TeV, 
(1) is the direct production, (2) - from Xcj decays, (3) - from ip' decays, (4) - sum of their all. For 
the different range in the rapidity: (a)- \y\ < 1.2, (b) - 1.2 < \y\ < 1.6, (c) - 1.6 < \y\ < 2.4 
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lfil, VS = 7TeV, 



FIG. 8: Prompt J/if) transverse momentum spectrum from LHCb Collaboration 
(1) is the direct production, (2) - from Xcj decays, (3) - from tp' decays, (4) - sum of their all. For 
the different range in the rapidity: (a)- 2.0 < \y\ < 2.5, (b) - 2.5 < |y| < 3.0, (c) - 3.0 < \y\ < 3.5 
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lfil, VS = 7TeV, 



FIG. 9: Prompt J/ip transverse momentum spectrum from LHCb Collaboration 
(1) is the direct production, (2) - from Xcj decays, (3) - from tp' decays, (4) - sum of their all. For 
the different range in the rapidity: (a)- 3.5 < \y\ < 4.0, (b) - 4.0 < |y| < 4.5 
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FIG. 10: The relative contributions of the color-singlet (curve 1) and color-octet (curve 2) pro- 
duction mechanisms to the prompt J/i/j transverse momentum spectrum at the \/~S = 7 TeV, 
1.5 < \y\ < 2.0. 
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